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In the changepoint problem, we determine when the distribution observed has changed to another 
one. We expand this problem to the quantum case where copies of an unknown pure state are 
being distributed. We study the fundamental case, which has only two candidates to choose. This 
problem is equal to identifying a given state with one of the two unknown states when multiple 
copies of the states are provided. In this paper, we assume that two candidate states are distributed 
independently and uniformly in the space of the whole pure states. The minimum of the averaged 
error probability is given and the optimal POVM is defined as to obtain it. Using this POVM, we 
also compute the error probability which depends on the inner product. These analytical results 
allow us to calculate the value in the asymptotic case, where this problem approaches to the usual 
discrimination problem. 

PACS numbers: 02.20.-a, 03.65.Ta, 03.65.Wj, 05.45.Tp 



I. INTRODUCTION 

The changepoint problem, which is studied in many 
fields, originally arose out of considerations of quality 
control. When a process is "in control," products are 
produced according to some rule. At the unknown point 
the process jumps "out of control" and ensuing products 
are produced according to another rule. It is necessary 
to determine the changepoint. 

In the classical case, we observe sequentially a discrete 
series of independent observations Xq,X±,... whose 
distribution possibly changes at an unknown point in 
time. It is assumed that independent random variables 
Xq, Ad, . . . , A„_i are each distributed according to a dis- 
tribution and the remaining independent random vari- 
ables X v , X„ + i, . . . are each distributed according to an- 
other distribution. Our purpose is to detect the change- 
point V Q. 

We extend this problem to the quantum setting where 
copies of an unknown pure state are being distributed 
in discrete time. Consider now the device distributing 
copies of an unknown pure state. This device has the 
unknown changepoint and distributes copies of another 
unknown pure state after the changepoint. Our goal is 
estimating the changepoint by observing all copies the 
device distributed. 

In this paper, we deal with the fundamental case, that 
is, we have only two candidates for the changepoint. We 
assume that the device distributes unknown pure states 
p t on the d-dimensional space for the discrete time t = 
0, 1, 2, . . . , t 3 . The state p t is changed at the changepont 
t c . That is, the states po, . . . , pt a -i are identical, and the 
other states pt c , . ■ . , pt 3 are also identical. Further, we 
assume that there are two candidates of the changepoint 
t\ and t%. 



p ® ■■■ ® p ® p. 



System 1 SystemO System2 

t = '=1, t - t, t = t. 
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FIG. 1: Systems 1, 0, and 2 



In order to analyze this problem, we introduce System 
1, System and System 2 to denote the systems that are 
the composite systems corresponding to the time period 
< t < t\, t\ < t < ti and t% < t < ts, respectively, as is 
explained in Fig. Q] Our task is then choosing the correct 
changepoint t c by using all three systems. This problem 
is equal to decide whether the state in System coincides 
with the state in System 1 or the state in System 2. 

We derive the optimal POVM and the minimum of the 
averaged error probability. This minimum probability is 
already obtained when the numbers of copies of states 
in System 1 and System 2 are the same. Sentis et al. 
studied this problem in the qubits case, i.e., in the case 
of d = 2 [2] and A. Hayashi et al. discussed it when 
System only has one copy Q. Our result concerns the 
general case that has no restriction for the numbers of 
copies and the dimension in systems. 

We also compute the error probability under the ap- 
plication of our optimal POVM, not averaged, which de- 
pends on the inner product of two states in System 1 and 
System 2. As the numbers of copies in three systems ap- 
proach infinity, this error probability clearly approaches 
exponentially unless the inner product is 1. Hence, the 
convergence speed can be measured by the exponential 
decreasing rate. The exponential decreasing rate seems 
related to the quantum Chernoff bound, which gives the 
optimal exponential decreasing rate of two state discrim- 
ination when a large number of copies of the unknown 
state are available [5] , while this relation has not been dis- 
cussed. In this paper, using the above analytical result, 
we calculate the exponential decreasing rate and clarify 
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the relationship with the quantum Chernoff bound. 

The paper is organized as follows. In the next Section 
we give the optimal strategy for the minimum averaged 
error probability. The optimal POVM is described by 
using the representation theory for easy calculation in the 
following sections. In Section III we obtain the minimum 
error probability represented as a function of the number 
of copies in each three systems and the dimension of the 
state space. Using the optimal POVM, in Section IV 
we compute the error probability which depends on the 
inner product of two states in System 1 and System 2. 
In Section V we consider the asymptotic behaviors of the 
minimum averaged error probability in several scenarios. 
In Section VI, we finally compute the convergence speed 
of the error probability. Some brief conclusions follow 
and we end up with a technical appendix. 



II. THE OPTIMAL POVM 

In order to treat the decision problem of the change- 
point, we denote the numbers t\, t 2 — ti, — t 2 + 1 by 
N\, M, N 2 , respectively. Then, we have Ni (N 2 ) copies 
of unknown pure states in System 1 (2). In the follow- 
ing, we denote the unknown state on System 1 (2) by p\ 
(p 2 ), which is a pure state on the d dimensional vector 
space C d . System has M copies of the unknown state 
p that is guaranteed to be either one of p\ and p 2 . Note 
that we assume that N\ < N 2 which loses no generality 
of this problem. Our purpose is to identify the state p 
with one of the two states by using all systems. This is 
equal to distinguishing two states, pf N 
and pf Nl ® pf A 
with equal probability. 

In this decision problem, we apply two-valued POVM 
{Ei,E 2 }, in which E\ (E 2 ) corresponds to the decision 
t c = ti (t c =t 2 ). Since E x = I - E 2 , our POVM can be 
described by a Hermitian matrix < E 2 < I, where / is 
the unit matrix. Then, the error probability is defined as 



pf M ® pf N * 
pf N2 , which are assumed to occur 



PM,Nx ,n 2 (Pi iP2,E 2 ) 
= i Tr[(pf N * ® pf M+N ')E 2 ] + ~ Tr[(pf M+JVl ® pf^E,] 

=±Tr[(pf N i ® pf M+N *)E 2 ] 

+ 1 -Tr[(pf M+N ^pf^)(I-E 2 )}. (I) 

Now, we assume that pi and p 2 are independently dis- 
tributed according to the unitary invariant distribution 
p® d on the set O^ of pure states on the d dimensional 
vector space C d . By using the POVM {/ - E 2 , E 2 }, we 
can define the averaged error probability as 



Pm,Ni,n 2 ( e z) 



l&d J@d 



PM,m , n 2 (Pi , P2 , E 2 )p Qd (dpi )pe d (dp 2 ) . 



(2) 



Here it is very helpful to use the following formula for the 
integral of the tensor product of L identically prepared 
pure states 



(3) 



where a is a pure state and II is the projector onto the 
totally symmetric subspace of (<C d )® L . 

By using this formula, the averaged error probability 
reads 

Pm,n u n 2 ( e 2) 
1 



f + Tr[( )E 2 ] 



where A\ and A 2 arc defined as follows: 



A, 



A 2 = 



N 1 +d-l\/M + N 2 +d-l 
d-1 J V d-l 



M + Nx + d-1 
d-l 



N 2 +d-l 
d-l 



(4) 



(5) 



(6) 



Note that A\ < A 2 , and the equation holds if and only 
if Ni = N 2 . 

Eq. (j4]) guarantees that the optimal strategy to mini- 
mize the averaged error probability is given by the Her- 
mitian matrix 

p _ ( In! ® ht+N 2 Im+n,. ® In 2 n1 
E m ,n 1 ,n 2 = { -. -. < 0}, (7) 



A, 



An 



where {A < 0} represents a projector onto the 
eigenspaces with negative eigenvalues of A. That is, plug- 
ging Eq. ([7]) into Eq. one obtains the minimum 
averaged error probability. 

In order to compute the minimum averaged error prob- 
ability, we deform the expression of the optimal POVM 
by using the tensor product representation of the unitary 
group U(d). 

Any irreducible representation of the unitary group 
\j(d) is characterized by a Young diagram A = 
[Ai, A 2 , . . . , Xd] and is denoted by U\. We use the short- 
hand notations X\ and [A 1? A 2 ] to denote [Ai, 0, 0, . . . , 0] 
and [Ai, A 2 , 0, 0, . . . , 0]. Note that Ul means the totally 
symmetric subspace of (C d )® L . The dimension of ,\ 2 ] 
is given as: 



dimWi 



[Ai,A 2 ] 



(Ai + d - 1)!(A 2 + d- 2)!(Ax - A 2 + I) 
(d-l)!(d-2)!(Ai + l)!A 2 ! 

(8) 

In our problem, the total system size is N = M + N\ + 
N 2 , and the total tensor product space (C d )® N can be 
decomposed to 



(C 



(9) 
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where Va corresponds to the multiplicity of the irre- 
ducible space U\. 

Since the tensor product space {$Z d }® N contains the 
two subspaces <8> Mm+n 2 an d Um+Ni ® Un 2 without 
multiplicity, these two subspaces have the form: 



U Nl ®Um+N 2 = ® k L U{N-k,k] ® C|u fc ), 



(10) 



u M+Nl ® u N2 = ®Z a o M+Nl ' N %N-k, k] ® C\v k ) (11) 

by using two normalized vectors: 

K)gV [JV - M] (0<k<Nt), (12) 
|«fc> G V_ fc , fc] (0 < fc < min(M + Aq, AT 2 )), (13) 

satisfying (u k \v k ) > 0. Since the dimension of Vat is one, 
the relation \uq) = \vq) holds. 

Letting I\ be the projector onto the space U\, one 
obtains from the equations (|T0l) and (fTTj) 

JVi 

JiVj ® Im+n 2 = J, I\N-k,k] ® |wfe)(u fe |, (14) 

fc=0 

min(M+ATi,JV 2 ) 

/m+ATi ® ijV 2 = X! V-fe.fe] ® KX u fe|- ( 15 ) 

Here, we note that the ranges of summations are different 
from each other. Using these equations, one has 

In ± <8> Im+n 2 Im+n-i ® In 2 



At 

^[jv-fe,fe] 



-4 2 

\u k )(uk\ \v k )(Vk\ 



At 



A 2 

\v k ){vk\ 



(16) 



min(M+Afi,AT 2 ) 

^[jv-fe,fe] 

fc=JVi+l 



Since |ufc) and |i;fc) are linearly independent in the range 
1 < fc < Nt (we show in Eq. (f28|)). there exists only one 
negative eigenvalue of \ Uk ^ Uk \ — l^*) v"» I _ Using the nor- 
malized eigenvector \wk) € Vrjv-fc.fc] with this eigenvalue, 
we therefore can write the optimal POVM as 



Ni 



E M ,N U N 2 = ^2l[N-k,k] ® kfcX^fcl 



fe=l 



min(M+Af 1 ,Af 2 ) 

X] I[N-k,k] ® |wfc>(«fc| 

fe=iVi+l 



(17) 



III. THE MINIMUM AVERAGED ERROR 
PROBABILITY 

In this section, we will compute the minimum averaged 
error probability. Plugging equations (TTBl and (TT7|) into 



Eq. Q as i?2 = Em,Ni,N31 one obtains 

PM,N t ,N 2 ( E M,Ni,N 2 ) 



iVi 



y^Tr[J[jv_ fc|fc ]](wfcl I 



|lifcX«fc| kfc)(Vjt| 



fe=i 



A 2 



^ mm(M+Ni,N 2 ) 
— £ Tr[V_ fc , fc] ] + l 

k=Nt+l 



\wk) 



(18) 



When arbitrary two real non-zero constants C\ and C 2 
and two linearly independent normalized vectors |a) and 
| b) are given, the unique negative eigenvalue of a ^ a 

\b)(b\ ■ ■ i 

J is given by 



Ci 



c 2 



C 2 - Ci - ^{C 2 -Ctf +4CiC 2 (l - |<a|6)| 2 ) 



(19) 



2CiC 2 

Therefore, the eigenvector \wk) associated with the neg- 
ative eigenvalue satisfies 



At 



An 



A 2 -At~ ^2-^)2+4^2(1-1(^1^)12) 



2AtA 2 

We also obtain the following equations: 

TM-fyv-fc.fc]] = dimW[jv_ feife ], 

fc=0 

min(M+A' 1 ,A r 2) 
fe=0 

Using these equations, we can write the minimum aver- 
aged error probability as 

Pm,n u n 2 (Em,Ni,n 2 ) = 



(20) 

(21) 
(22) 

(23) 



At + A 2 



An 



E ^ U A N . kM V(A 2 - AtY + AMM1 - \{u k \v k )V) 

fc=0 1 2 

(24) 

Our remained task is to calculate the inner product 
(u k \v k ). When we denote the highest weight vector of 
the space ZY[jv-fe,fc] by \[N— fc, fc] d ), (u k \v k ) is equal to the 
inner product of | [N - k. k] d ) \u k ) and \[N - fc, k] d )\u k ). 
We can assume d — 2 without loss of generality since the 
inner product does not depend on the dimension. Let us 
fix some notations as follows: 

_ M _ Nt _ N 2 

M + Nt M + N 2 _ N , 
M01 = 7; ,^02 = 7: ,fi= "7; fc- 



(25) 
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Using Wigner's 6j-function Q, we then can write 

(Uk\Vk) 

= (-ir+^ 2+ V(2Moi + l)(2/i 02 + 1) ^ t "° 1 } 

f- 1 P02J 



(26) 



Moreover, Wigner's 6j-function can be computed as 



Mi Mo A*oi 
M2 M M02 



(27) 



V(2moi + 1)(2a*02 + 1)^ ( M t Wl )( M r 2 ) 
(Appendix 1). Thus, one obtains 



(ufc|ufc) = 



\ 



(28) 



and in order to denote this value we use the notation (f>k 
satisfying 



cos <pk 



ft) ft) 



\ (T)(T) 



\Uk\Vk) 



Therefore, the minimum averaged error probability can 
be written as 



Pm,n 1 ,n 2 { e m,n u n 2 



Ao 



y-, dimZY[jv-fc ,fe] 



fc=0 



^1^2 



(A 2 - Ai) 2 +4AiA 2 sin 2 fe 



In the case of d = 2, since 

i4 1 = (JVi + l)(M + J\T 2 + l) 

A 2 = (M + JVi + l)(JV 2 + l) 

A 2 - Aj = M(JV 2 - iVi) 

dimW[jv_ fe)fc] = M + Ni + N 2 - 2k + 1, 



(29) (|30l) is calculated to the following way. 



(30) 



Pm,n u n 2 ( e m,n u n 2 ) 



(JVi + l)(M + JV a + l) 



(Af + JVi + 1)(JV 2 + 1) f-' (JVi + 1 M + N 2 + 1 (M + JV X + l)(JV a + 1) 



M + Ni + N 2 - 2k + 1 



'M 2 (iV 2 - TVx) 2 + 4(JV X + 1)(M + iV 2 + 1)(M + JVi + l)(iV 2 + 1) 



1 - 



ATx^M + iVi - fc)! \ / N 2 \{M + N 2 ~k)\ 



(M + Ni)\(N 1 - fc)! J \ {M + N 2 )\(N 2 - k)\ 

(31) 



When iVi — N 2l the equation A% — A 2 holds and this 
probability is concretely computed as 



PM,N x ,N x {EM,N lt Nx) 

l\ {d - l)Nx\{M + Nt)\ 

~ 2 ~ (Ni+d - 1)!(M + N 1 + d-l)\ 

x + 2Ni -2k + l) 

k=0 

(M + 2N 1 -k + d- l)\(k + d-2)\ 
X (M + 2Ni -k + l)\k\ 



x Wl- 



N L \{M + Nx - k)\ 
{M + N 1 )\(N 1 - k)\ 



(32) 



Moreover, plugging d = 2 into this equation, we have 



fe=0 



M + 2N t - 2k + 1 
(N 1 + 1)(M + N 1 + 1) (33) 



x A /l-( 



N^iM + Ni-ky. 
(M + iVi)!(7Vi - fc)! 



This result coincides with the result by Sentis et 



IV. THE ERROR PROBABILITY WITH THE 
OPTIMAL POVM 

In the previous section, we obtained the averaged er- 
ror probability when the two candidate states are dis- 



tributed independently. However, the unknown states 
Pi and p 2 do not necessarily obey the uniform distribu- 
tion. In order to treat the performance of our Optimal 
POVM in a more general setting, we consider the er- 
ror probability for the given two pure states p\ and p 2 
when the optimal POVM is applied. This error proba- 
bility depends on the inner product of two candidates, 
i.e., q = Tr[pip2]- In the following, we calculate the error 
probability given by Eq. ([1]) in the case of the optimal 

POVM {I — Em,Ni,N-i,Em,N-l,Nz}- 

Theorem 1 The error probability with the optimal 
POVM {I — Em,n 1 ,N2^m,n 1 ,n 2 } can be written as 

PM,Nx ,N 2 (Pi 5 92 , Em,N! ,n 2 ) 



1 N1 



k=0 



Pk + Qk 



(A 2 - A x )( P k -Q fc ) + 2sin 2 <t> k (A x P k + A 2 Q k ) 
(A 2 - Ai) 2 +AAxA 2 sin 2 4> k 



(34) 



where P k and Q k is given as follows: 
(V-2fc + l)Vi!(M + N 2 )\ 



Pk 



Ni-k 

x E 

1=0 



(N-k + l)\k\ 
N 1 -k\fM + N 2 -k + l\ , X| . , 



(35) 



(N - 2k + 1)(M + iVi)!JV 2 ! 



(N-k + l)\k\ 



N 2 -k 



E { M + N 1 - k + iYN 2 -ky {i _ q)N2 _ l 



i 



Here we have defined 0=1. 
When N± = N 2 , we can write 

Pm,n 1 ,n 1 (pi, Pi, Em.Nx^Ni) 



I 



Ni 



k=0 



(M + 2Ni ~2k + l)Vi!(M + Ni)\ 
(M + 2Vi-A: + l)!fc! 



xyi- 

Ni-k 



( N 1 l{M + N 1 - k)\ 



x 



\(M + Ni)\(N! -k)\ 



1=0 



Ni-l 



(37) 



Proof. Let us start by defining the notation P k and 
Qk as 



P k = Tr[I [N _ kM ® \u k )(u k \(pf ^ ® p® M+N2 )] 
(0<k< Ni), 



(38) 



} k = Tr[I [N _ kyk] (g) \v k )(v k \(pf 



M+Ni „ ®N 2 



(39) 



(40) 



(0 < k < mm(M + N 1: N 2 )). 
We note that the followin equation holds: 

Ni mm(M+Ni,N2) 
k=0 k=0 

Since an arbitrary pure state a satisfies 
one obtains 

Tr[( P ri®pf M+N *)E M!NuN2 ] 
= Tr[(I Nl ® Im+nM Nl ® pf M+N2 )(I Nl ® /m+jv 2 ) 

(42) 

Plugging equations (fT4|) and (|17[) into Eq. (j42|) . one has 

Ni 

Tr[(p? Wl ® ^^M^.ivJ = £ IKK>| 2 ^- 

fc=i 

(43) 

In the same way, using equations (|15|) and (|42j) . we obtain 

ir l(Pl J-E^MjNi ,N 2 \ 

Ni mia(M+Ni,N 2 ) 

= ^\{vk\w k )\ 2 Q k + Yl ® k - 

k=l k=Ni+l 

When two arbitrary normalized and linearly inde- 
pendent vectors | a), 1 6) and two positive real numbers 
Ci,C 2 > are given, the normalized eigenvector |— ) with 
the unique negative eigenvalue of a ^ a — MiM satisfies 



(36) the following equations: 



Ci c 2 



\(a\-)\ 2 = l 



1 - 



C 2 -C 1 + 2C 1 (l-\(a\b)\ 2 ) 
V(^2-C 1 ) 2 + 4C 1 C 2 (l-|(a|6)|2) 

(45) 

C 2 -C 1 ~2C 2 (l-|(a|b)| 2 ) 
V(^-C Y 1 ) 2 + 4C 1 C 2 (l-|(a|6)|2) 

(46) 



Applying equations (|45p and (|4l)|) to the case of \a) — 
\u k ), \b) = \v k ), we have 



\(u k \w k )\ 2 = -(1 



\(v k \w k )\ 2 = -(I 



A 2 


-^4 


2A\ sin 2 (j) k 


y/(A 2 - 


-^i) 2 


+- 4AiA 2 sin 2 <p k 


A 2 


-A 1 - 


2A 2 sin 2 fc 


J(A 2 - 


-A^ 


■f 4AiA 2 sin 2 (j) k 



), 

(47) 
)■ 

(48) 
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for 1 < k < Ni, 

Using these equations, we can write the error proba- 
bility as 

PM,Ni ,n 2 (pi i Pi i Em.Nx ,n 2 ) 



We can calculate the Clebsch-Gordan coefficients @ as 
(N-2k + l)(Ni-k)l(Ni-l)\ 



n ®N 2 \j? 1 
P2 )- tj M,N 1 ,N 2 \ 



(N - k + l)!(iVi - k - 1)1 

(M + N 2 y.(M + N 2 - k + iy. 
X {M + N 2 - k)\k\l\ ' 



(54) 



1 Nl 



JVi Ni 

1 +e \(u k \wk)\ 2 p k - E i(^K>i 2 Qfc 

k=l k=l 
miniM+Nt.N?) 

E ^ 

Pfc + Qfc 



Denoting pi = |0i)(0i|, we obtain \{<pi\ t)| 2 = 3- 
Thus, 

+0 



Aq 



|(<H t)\ 2l \(^i\i)\ 2{Nl - l) 



(55) 



I 



(A 2 - Ai)(P fc - Q fc ) + 2 sin 2 M^iPk + MQk) 



(A 2 - AiY+AAxAz sin 2 



(49) 



P 



Therefore, we can write Pk as 
(N - 2k + l)Aq!(M + N 2 )\ 



Now we turn our attention to computing Pk- We can 
assume d — 2 since Pk does not depend on the dimension. 
By using Clebsch-Gordan coefficients and the notations 
given by (|25|) . the projector in Eq. (|38|) can be written 
as 

I[N-k,k] ® \Uk)(Uk\ 
N-lk Ni M+N 2 

= E E E \(p '■ p - l \pi : pi - ^ P02 ■■ P02 - j)\ 2 

1=0 i=0 j=0 

x \fii : p ± - i)(px : p x - i\ <X> \p i : p Q i -j)(p 02 : p 02 - j\ 

(50) 

In the following, we fix the notation | ( | \) ) to 
denote the vector in the space C 2 whose weight is \ (— |). 
Without loss of generality, we can assume that p 2 is J t 
) (t I • Then, we can write 

pf M+N 2 = |Mo2 . Mo2)(Mo2 . Mo2| = ! t)(t ^M+N^ (51) 

Plugging equations (f50|) and (|5T|) into Eq. (f38|) , we have 



E 

2=0 



(N — k + l)!fc! 

iVi - fc\ /M + N 2 - k + 1\ j 



g'Ci-g)^" 1 . 



(56) 



where we have defined 0° as 1. 
In the same way, one obtains 



Qi, 



(N — 2k + 1)(M + Aq)!iV 2 ! 



(N -k + l)\kl 



N 2 -k 



Nt-k 



p k = e <w : mi - * - iipf^iw pi — k—i) 



1=0 



(52) 



x |(/x : /x - Z : /xi - fc - Z;/x 2 : M02) 



Moreover, converting the variable Z into N% — k — I, this 
can be written as 



iVi-Jfc 



,Aq iVi 



2=0 



xE( M+Ar r" +z )rr>' (i -^- 

(57) 
□ 



V. LIMIT OF THE MINIMUM AVERAGED 
ERROR PROBABILITY 

When the both numbers of copies in System 1 and 
System 2 approach infinitely large, we have perfect 
knowledge to determine the states in the two sys- 
tems. In this limit, by using Eq. (I32[) . the probability 
Pm,n 1 ,n 1 { e m,n 1 .n 1 ) can be written as 

lim Pm,Ni,Ni( E M,Ni,Ni) 
iVi— >oo 

_ 1 

~ 2 



1 - 2(d- 1) / xy/l-x 2M {l~x 2 ) d - 2 dx 



..N N AT JVi xVi , . M + N 2 M + N 2s 
x\( Y -k: Y -N 1 + l\—:-— + l;^^:^^ l 

(53) 



I 2 - =2 



l-(d-l) / Vl-x M (l~x) d - 2 dx 
Jo 



(58) 
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where we have defined x = ^ and used the Euler- 
McLaurin summation formula. The case of d = 2 coin- 
cides with (18) in Santis et al 0, and the case of M = 1 
coincides with (41) in A. Hayashi et al Q. 

This result could be easily anticipated from the mini- 
mum error probability of the discrimination problem 0] ■ 
Recall that the minimum error probability given M iden- 



tical copies is 



1 



Assuming that 



(Tr[pip 2 ]) M 

pi and p 2 are distributed according to p,Q d independently, 
the average is given by 



e d Je d 
1 

" 2 
1 

" 2 



1 - Jl - (Tr[p lP2 }) 



M 



Pe d {dpi)pe d (dp 2 ) 



1 — 2(d — 1) / sj\ - (cos(9) 2M (sin 
o 
i 



\2d-3 



l-(d-l) / y/l-x M (l-x) d - 2 dx 

Jo 



cos OdO 



(59) 



Therefore, our optimal measurement can achieve the av- 
erage performance of two state discrimination under the 
limit N\ = N-2 — > oo. 

Next, we turn our attention to the complementary 
case, that is, the number of copies in System is in- 
finitely large. By using Eq. (|30[) . the minimum averaged 
error probability in this limit can be computed as 



Since this value is the same as the expression in Eq. (|60l) , 



the POVM {E 1 = I x ® (I 2 
realizes the optimal performance. 



),E 2 = h 



} 



VI. EXPONENTIAL DECREASING RATE OF 
THE ERROR PROBABILITY 

When every number of copies N±, N 2 and M is in- 
finitely large, that is, the states in three systems are per- 
fectly known, the error probability approaches zero unless 
q = 1. In order to treat the convergence speed, we focus 
on the exponential decreasing rate of the error probabil- 
ity when the optimal POVM Em,n 1 ,n 1 is applied. For 
simplicity, we assume that the numbers Nx, N 2 of copies 
in Systems 1, 2 increase in proportion to the number M 
of copies in System 0. When the proportional constant 
is given to be a > 0, we have N\ — N 2 — aM. Thus, 
using the real numbers 

_ (M + 2aM -2k+ l)(aM)\{M + aM)\ 



2(M + 2aM - k + l)!fc! 



1-1 1 



(aM)\(M + aM — k)\ 
(M + aM)\(aM — k)\ 



(63) 



M 



1 



(60) 



Thus, the average is computed as follows: 

/ \ (Trf/vD^ 2 Pe d {dp)pe d (dp) 
e d J&d 1 

= (d-l) [ 2 (cose) 2N2 ( S m9) 2d ' 3 cos 



x N2 {l-x) d -'dx 



(61) 



d-1 f 1 



(62) 



1 



D k i = 



aM - k\ (M + aM - k + l\ ^ _ 



Note that this result is independent of N\. 

In this limit, we have perfect knowledge of the pure 
state p in System and this problem is equal to distin- 
guishing two states p® Nl ® p® N2 and p® Nl ® p® N2 in 
the composite system 12. This problem can be regarded 
as a generalization of state comparison Q. As is shown 
in the following, the minimum error probability for these 
two states can be obtained with a POVM whose elements 
are {E x = h®(I 2 - p® N2 ), E 2 = h®p® N2 }, where h (I 2 ) 
is the unit matrix on (C^)®^ 1 {{C d )® N2 ). Here, £i(£ 2 ) 
corresponds to the guess p® Nl <g> p® N2 (p® Nl ® p® N2 ). We 
then can write the error probability as 

\ ^[p® N ^p® N2 E 2 ]+ l - TryP^top 9 "*^] = \ (Tr[pp]f 2 , 



I J\ I 
we can write from Eq. (|37p 

Pm 



(64) 



aM aM—k 



,aM,aM(Pl, Pi, %,aitf,atf) — ^ ^ CfcDjy. (65) 



k=0 1=0 



which 



The convergence speed is represented as 

nm ~TT l°&PM,aM,aM(pi, P2, EM,aM,aM), 

M— ¥00 M 

can be deformed as 



Um T7 \°§>PM,aM.aM(Pl, P2, Em.oM ,olm) 
M-s-oo M 



aM aM—k 



k=0 1=0 



M-s-oo M 

um TT lo S( max max C k D k i). 



(66) 



Moreover, using the approximation formula \/l — x 
1 — hx( when x <C l)and the Stirling approximation n\ 
n n e~ n \ // 2~7m, one obtains 



lim — \0gPM.aM,aM{pl,P2 1 EM,aM.aM) 

M-s-oo M 

= min min /i(/3, 7), 

0</3<a0<j<a-l3 
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where we have defined as (3 



M ' ' 



/ l (/3, 7 ) = (a-^- 7 )log(a-^-7) 
- (1 + a - + 7) log(l + a - /3 + 7) + (a - log(a - 0) 
-(1 + a-fi) log(l + a - /3) - 3a log a 
+ (1 + a) log(l + a) + /3\og/3 
+ (1 + 2a - p) log(l + 2a - /3) 
+ 2 7 log 7 -7logg- (a-7)log(l - g). 

(68) 

There is the unique root of ^4 = = in the range 
0</3<a, 0<7<a — /? and we use (/3i, 71) to denote 
it. These can be calculated as 



71 



q(a 



1 



A = 



(2a + 1) - ^(20 + 1)2-4(02 -a 7l ) 



(69) 



(70) 



One can agree that h(J3i, 71) is the minimum of the func- 
tion h in the range < (3 < a, < 7 < a — (3 due to the 
following equations: 

d 2 h, , 
iH(A7)>0, 



lim 



9/3 2 
<9 2 /i. 
97 2 

' , ^, 7 )=lim^ l7 ) 



v2 "(^7)>0, 



(71) 



dp 



+0 97 



dh, n \ ,. 9h. n 



When a is sufficient large, we can write 

1-9 



MA, 7i) =-logg- 



q(a-l) + 2 



(72) 



In fact, as is numerically demonstrated in Figs. [5] and 02 
— log q — q (a-i)+2 weu approximates ft.(/3i, 71) when a is 
large. Therefore, we obtain the convergence speed of the 
error probability, 



um ~TT logPM.aM.aAf (Pl> P2, ^M.aM.aAf) 
Ai— foo 



= -logTr[p 1 p 2 ] - 



1 - Tr[pi,p 2 ] 
Tr[p 1 p 2 ](a-l)+2 



(73) 



In the discrimination problem of two pure states [J], 
when the number of copies of the state to be identified is 
infinitely large, the convergence speed is given by 



lim — loi 

M-»oo M 



1-Jl-fTtM) 



M 



-logTr[/9ip 2 ]- 
_ (74) 

This is called the quantum Chernoff bound [5[ and equal 
to the limit of Eq. (|73|) as a — > 00. This fact means 
that the performance of our optimal POVM is close to 
that of the optimal POVM in the sense of quantum state 
discrimination. 



Convergence Speed 
3.0 F 
2.5 
2.0 
1.5 
1.0 
0.5 



FIG. 2: (normal line) ft(/8i, 7 i), (dashed line) —log 
q (a-i)+2 ' ( thick line ) ~ [o S1 for a = 5 



Convergence Speed 
0.7 e 



10 12 14 



FIG. 3: (normal line) h(fii,"/i), (dashed line) — logg 
q (a-i)+2 ' ( thick line ) - log g for g = 0.5 



VII. CONCLUSIONS 



We have studied the changepoint problem in the quan- 
tum setting, where our task is to choose the correct 
changepoint between two candidates. This problem is 
equal to discriminating two unknown general states when 
multiple copies of the state are provided. We have ob- 
tained the minimum averaged error probability, Eq. (|30p . 
Our result of special cases coincides with the results by 
and ||. However, when arbitrary numbers of copies 
of general pure states are given, we have calculated it 
for the first time. Moreover, we have first calculated the 
non-averaged error probability Eq. (|34p . This depends 
on the inner product. As could be anticipated, when 
the numbers of copies of candidates are infinitely large, 
we recover the average of the usual discrimination prob- 
lem. We have also paid attention to the exponential 
decreasing rate and shown the convergence rate of the 
non-averaged error probability approaches the quantum 
Chernoff bound. 



9 



Acknowledgment 



This is deformed as 



The authors thank Professor Leong-Chuan Kwek for 
pointing out the importance of the changepoint problem. 
They are grateful for Dr. Gen Kimura and Mr. Wataru 
Kumagai to helpful discussion concerning qunautm infor- 
mation. The authors are partially supported by a Grant- 
in- Aid for Scientific Research in the Priority Area "Deep- 
ening and Expansion of Statistical Mechanical Informat- 
ics (DEX-SMI)," No. 18079014 and a MEXT Grant-in- 
Aid for Young Scientists (A) No. 20686026. The Centre 
for Quantum Technologies is funded by the Singapore 
Ministry of Education and the National Research Foun- 
dation as part of the Research Centres of Excellence pro- 
gramme. 



Appendix A: Calculation of Wigner's 6j-function 

Let us consider the calculation of the Wigners 6j- 

function 1 a , ^ C r \. Let us define some notations as 
Id e fj 

ai = a+ b + d + e,a 2 = a + c + d+f,a 3 = b + c+ e + f, 
Pi = a + b + c, fi 2 = a + e + f, 
P 3 = b + d + /, #4 = c + d + e, 

(Al) 

and let us define Ai,A 2 and A 3 to be the smallest, mid- 
dle, and largest values of ai,«2 and a 3 and B\,B 2 ,B 3 
and B4 to be the smallest, second smallest, second 
largest, and largest values of @i,P 2 ,P 3 ancl Pa- When 
£? 4 = Ai, from the formula in p| we can calculate as 



a b c 
d e / 



(-1) S4 
1 



-B4 — B3 

n 



[B 3 + 1 + i){A 3 - Ai + i) 
Ll {A 3 -Bi + i)(A 2 -B 2 +i) 



(Bi + 1)(B 2 + 1) 

a 2 -a 



n 



{A 1 - B l + i)(Bj ~B 2 + i)(Bj -B 3 +i) 
(Ai +B X - A 2 + i)(Ai +B 2 -A 2 + i)i 



(A2) 

We now compute the Wignerfs 6-j funcion 

mi w> Moil inE 
M2 m M02 J 



In case < k < M 



Ml Mo Moil 
M2 M M02 J 

' V( 2 ^oi + l)(2Aio2 + l)\ 



(A5) 



In case M + 1 < fc < Ni 
In this case, we can write as 

Ai=N- k, A 2 = M + N — k, A 3 = N 

B x = M + N U B 2 = M + N 2 ,B 3 = N -k,B 4 = N - k. 

(A6) 

Plugging these into Eq. (A2), one obtains 

= {-l'f- k 



Mi Mo M01 
M2 M M02 



(M + Ni+1)(M + N 2 + 1) 



M 

n 



(Ni -k + i){N 2 -k + i) 
(Ni+i)(N 2 + i) 



Since M + 1 < k, one has 



Thus, 



M 

n 

i=l 



(Aq - fc + i)(JV 2 - fc + i) 
(Aq+i)(A 2 + z) 



2 


[Mi 


Mo 


Moil 




iM2 


M 


M02 J 



v(2moi+i)(2mo2+i)^ rr i )( M r 2 ) 

also holds. 



(A7) 



(A8) 



(A9) 



In this case, we can write as 
Ai = N - k, A 2 = N, A 3 = AI + N - k 
B 1 = M + N U B 2 =M + N 2 ,B 3 = N -k,B A = N - k. 

(A3) 

Plugging these into Eq. (A2), one obtains 



Mi Mo M01 
M2 M M02 



(-1) 



N-k 



1 



n 



(M + Ni + 1)(M + N 2 + 1) 
{Ni-k + i)(N 2 -k + i) 



L (M + Ni-k + i){M + N 2 -k + i) 



(A4) 
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